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Abstract 

The Modified Chaplygiii Gas (MCG) model belongs to the class of a unified models of 
dark energy (DE) and dark matter (DM). It is characterized by an equation of state (EoS) 
Pc = Bp — A/p", where the case B = corresponds to the Generalized Chaplygin Gas (GCG) 



O ' model. Using a perturbative analysis and power spectrum observational data we show that 

.4—) , the MCG model is not a sucessful candidate for the cosmic medium unless B = 0. In this 

^ ' case, it reduces to the usual GCG model. 

PACS: 98.80.-k, 04.62.+V 

1 Introduction 

f""- ■ The cross of different observational results at cosmological level indicates that besides the usual 



expected contents of the cosmic budget, like baryons and radiation, there is a dark sector, 
with two components, dark matter and dark energy. In principle, dark matter is present in local 
structures like galaxies and cluster of galaxies, suffering consequently the process of gravitational 
collapse. In this sense, dark matter behaves much as like ordinary matter. However, it does 
emit any kind of eletromagnetic radiation. Dark energy, on the other hand, seems to remain 
C^ ' a smooth, not clustered, component, driving the accelerated expansion of the universe. This 

property requires a negative pressure. Many different models have been evoked to describe this 
dark sector of the energy content of the universe, going from the inclusion of exotic components in 
the context of general relativity theory to modifications of the gravitational theory itself, passing 
by other possibilities as the breakdown of the homogeneity condition. For a recent review, see 
reference [T]. 

A very appealing proposal to describe the dark sector are the so-called unified models. The 
protype of such model is the Chaplygin gas O [31 |4]. In the unified model dark matter and 
dark energy are described by a single fluid, which behaves as ordinary matter in the past, 
and as a cosmological constant term in the future. In this sense, it interpolates the different 
periods of evolution of the universe, including the present stage of accelerated expansion. The 
Chaplygin gas model leads to very good results when confronted with the obsevational data of 
supernova type la [5]. Concerning the matter power spectrum data, the statistic analysis leads to 
results competitive with the ACDM model, but the unified (called quartessence) scenario must 
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be imposed from the begining [H [7] . This means that the only pressureless component admitted 
is the usual baryonic one, otherwise there is a conflict between the contraints obtained from the 
matter power spectrum and the supernova tests. 

Many variations of the Chaplygin gas model have been proposed in the literature. One of 
them is the Modified Chaplygin gas (MCG) Model. The equation of state of the MCGM is 

Pc = Bp-Ap-^, (1) 

where B, A and a are constants. When i? = we recover the Generalized Chaplygin Gas (GCG) 
Model, and if in addition a = 1 we have the original Chaplygin gas model. The dynamics of 
the MCG model has been studied in the reference |8j, while a dynamical system analysis has 
been made in reference [9]. The evolution of the temperature function has been considered in 
reference [10]. Some background constraints were established in references |11| and |12| . The 
analysis of the spherical collapse was made in reference |13) , while a perturbative study, looking 
for some general features of the model, was carried out in reference [14]. In all these studies the 
viability of the model was concluded, but no one of them has exploited the observational data 
concerning the perturbative behaviour of the model. 

Our intention here is to test the MCG model against the power spectrum observational data. 
For the background tests, as those analyzed in references [11^ 112) . the MCG reveals to lead 
to competitive scenarios compared with the ACDM and the GCG models (to give just some 
examples). However, the constraints coming from the power spectrum data are, in general, 
much more crucial since it tests not only the background framework, but also the perturbative 
behaviour of the model. Many general constraints can be established on the parameters of the 
EoS ([1]) even not considering perturbations. For example, in the past the equation of state 
([1]) implies, when a is positive (a requirement necessary in a perturbative analysis in order to 
preserve a positive sound speed) that 

/'^(« - 0) = ^^- (2) 

In order not to spoil the usual primordial scenario of the standard model (in special nucleosyn- 
thesis), B must be smaller than 1/3, which may include negative values. On the other hand, as 
we will show below, the requirement that the sound speed of the MCG must be positive implies 
essentially that i? > 0. Hence, the admissible values of the parameter B seems to be around 
< i? < 1/3. These considerations will be strengthened through the power spectrum analysis 
to be made later in this paper: the matter power spectrum data can be fitted only if |i?| < 10~^. 
Hence, essentially the only configuration possible is that corresponding to the generalized Chap- 
lygin gas, with perhaps some possible very small deviations from it. In this sense, we can consider 
that the MCG model is ruled out when confronted with the power spectrum observational data. 
In next set we set out the general equations of the MCG model at background and perturba- 
tive levels. In section III we perform a numerical analysis comparing the theoretical results with 
the matter power spectrum observational data. In section IV we present our conclusions. 

2 Basic set of equations 

Our starting point are Einstein's equations coupled to a pressureless fluid, radiation and to the 
MCG fluid. They read, 

R^.. = 87rG|T--ig^,r-|+87rG|T;,-ig^,r'-|+87rG|r^^-i5M.r'=|, 

r-pt^V _ n. rp/lU _ rv rjiflU _ n. 

in ;^ ^ J C -jfj, ") r;/i " 

The superscripts (subscripts) m, r and c stand for "matter", "radiation" and "Chaplygin". We 
assume a perfect fluid structure for the cosmic medium as a whole and also for each of the 



components, 

^r = {pA + PA)u^A<i - PAg'"" , A = m,c,r . (3) 

Using now the flat Friedniann-Robertson- Walker metric (as suggested by the Seven-year WMAP 
data [15]), 

ds'^ = dt^ - a{tf[dx^ + dy^ + dz^] , 

and identifying all the background 4-velocities, Einstein's equations reduce to 
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In the above set of equations we have defined Ag 
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The perturbed equations in the synchronous coordinate condition can be established following 
closely the computation shown in reference [7] . We introduce fluctuations around the background 
quantities, 5^,^ = 5^,^ + h^^, p = p + 5p, p = p + 5p, u^ = u^ + 5u^. The bars indicate the 
background quantities. The synchronous coodinate condition implies /i^o = and 6u^ = 0. The 
final perturbed equations read (see also |16)). 
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where p and p stand for the total matter and pressure, respectively, 9 = 5u^^ and h = hk]^/a?. 
In terms of the components, we end up with the following equations: 
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with 6m = ^Um,i^ ^c = ^U'ci ^^^ Or = 5u\.^. 
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and remembering that Pm = ^Pm = 0, leading also to 0^ = up to an irrelevant non- 
homogeneous term, the set equations becomes 
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and /cq = 3000 hMpc. A Fourier decomposition of the spatial dependence of the perturbed 
quantities has been performed, k being the corresponding Fourier mode. 



3 Numerical analysis 

The matter power spectrum is defined by 



V 



(30) 



where 6k is the Fourier transform of the dimensionless density contrast 6m ■ In what follows we will 
use the matter power spectrum data of the 2dFGRS observational mapping |17j . We use the data 
in the range O.OlMpc^^h < k < 0.185Mpc~^h (here h is connected with the Hubble parameter 
today being defined by Hq = 100 hkm/Mpc.s, with no relation to the previous definition of 
the metric fluctuation) which corresponds to the linear regime which is being considered here. 
Hence, the covariance matrix is diagonal. 

One important aspect of the numerical analysis is to fix the initial conditions. In order to 
do so, we use a scale invariant primordial spectrum, with the BBKS transfer function |18) . The 
initial conditions are fixed following in general lines the prescription described in reference [19|. 
Even if a full statistical analysis may be performed using, for example, a bayesian statistical, 
such complete analysis is not essential to obtain the main result of this work. 
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Figure 1: In each panel we have fixed one of the parameters (a,As,S) and ploted the contours for which the 
speed of sound is equal to zero. Above (below) each dashed line we have u^ > (Va < 0) for differents values of 
a — (As) — (B) in the left- (center)- (right) panel. 

In any perturbative analysis, one crucial quantity is the squared sound speed. Positive values 
are assured if B and a are positive, see equation (21). It is possible to have also positive squared 
speed in the case a is negative, but in a very small range. The possibility to have v'^ > and 
S < is almost excluded, as it can be seen in figure 1. However, in performing the numerical 
computation we will consider this possibility. 

In figure 2 we plot the linear matter power spectrum comparing with the observational data 
for As = 0.95 and a = 10 (left pannel) and a = 1 (center pannel), with different values of B. 
In the right pannel the two-dimensional Probability Distribution Function is displayed, showing 
the contours at la, 2a and 2>a confidence levels. When a = 10, the case i? = is essentially the 
only one that fits the observational data, with a total x^ of about 17 (the statistical parameter 
that measure, to say in general lines, the quality of the fitting of the data by a theoretical model, 
the smaller the value of x^ the better the fitting). As it can be verified, when |-B| > 10~^ there 
is already a huge discrepancy between the theoretical results and the observational data. At 
the center pannel of the same figure, the only changing is a = 1. Now the better fitting is 
achieved by i3 = —2.7 x lO"'^, but with x^ ~ 37: the fitting is much worse than in preceding 
case. An extensive inspection of the different possibilities shows that a reasonable agreement, at 
la confidence level, can be find only around \B\ < 10^^, with a » 1 or a ~ 0. Such fine tuning 
in the EoS parameter ([T]) implies that the only "natural"' value would be i? = 0. This reduces 
the model to the the GCG model. Hence, the confrontation with the matter power spectrum 
data seems to rule out the MCG model. 



4 Conclusions 

The Modified Chaplygin Gas (MCG) model is an extension of the Generalized Chaplygin Gas 
(GCG) model, with the addition of a new term in the usual GCG equation of state which is 
proportional to the density. In this sense, it can be seen as a combination of the CGM model 
(where p = —A/p'^) and the so-called XCDM model (where p = Bp). In the present paper we 
have confronted the MCG model against the power spectrum observational data. In general, the 
MCG leads to good behaviour concerning the background tests. However, the power spectrum 
analysis implies to perform a perturbative study of the model, testing in this sense the deep 
content of the theoretical framework. In this perturbative analysis, a crucial quantity is the 
squared sound speed which must be positive in order to avoid instabilities. In the MCG model. 
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Figure 2: At left, power spectrum for the MCG model fixing As = 0.95 and a = 10. From top to bottom, the 
curves are MCG models with B = 10"*, B = 10"^, B = 0, B = -lO"'^ and B = -lO""*, respectively At center, 
power spectrum for a — 1.00, and for the same values as before of the parameters As and B. In the right pannel 
it is displayed the la, 2a and 3a confidence levels contours for the two-dimensional Probability Distribution 
Function for B and a. 



the positivity of the squared sound speed requires B and a to be also positive, even if small 
range of negative values of B and a can in principle be admitted. 

However, the confrontation of the MCG model with the power spectrum observational data 
rules out any significant departures from B = 0. In fact, only values such that \B\ < 10~^the 
model can fit the observational data. This result essentially rule out the MCG model, reducing 
it to the GCG model. Crossing these results with those of reference [SlI^, where power spectrum 
constraints on the GCG model were establisheds, only the quartessence version of the GCG 
model survives. In the quartessence formulation of the GCG model, the matter content is given 
only by the baryonic component. If this quantity is left free, we find a prediction of a universe 
dominated only by matter, in contradiction to the supernova analysis [5]. 

It must be remarked however that the conclusions obtained in the present work are restricted 
to a hydrodynamical representation of the dark energy component. If this component is alter- 
natively represented by a self-interacting scalar field, leading to the same background relations, 
the restrictions arising from the positivity of the squared sound speed disappears. If the MCG 
model can survive the observational tests in this alternative formulation is a question we intend 
to address in a future work. The main difficult comes from the fact that the usual GCG gas is 
connected a kind of Born-Infeld , which leads to essentially the same expression for the sound 
speed as in the hydrodynamical formulation [20]. Hence the restrictions found here for the hy- 
drodynamical formulation may remain. If such problem can be circumvented in the case MCG 
model is an open question. 
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